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Abstract. We conducted a 4 epoch 3 month VLBA proper motion study of H2O masers toward an intermediate- 
■£3 ■ mass class source i20 MMS1 (d =700 pc). The region of i20 contains at least 3 bright young stellar objects 

at millimeter to submillimeter wavelengths and shows three pairs of CO outflow lobes: the brightest source 
MMS1, which shows an extremely high velocity (EHV) wing emission, is believed to drive the outflow(s). From 
milli-arcsecond (mas) resolution VLBA images, we found two groups of H2O maser spots at the center of the 
K^j ' submillimeter core of MMS1. One group consists of more than ~ 50 intense maser spots; the other group consisting 

of several weaker maser spots is located at 18 AU south-west of the intense group. Distribution of the maser spots 
in the intense group shows an arc-shaped structure which includes the maser spots that showed a clear velocity 
gradient. The spatial and velocity structures of the maser spots in the arc-shape did not significantly change 
through the 4 epochs. Furthermore, we found a relative proper motion between the two groups. Their projected 
separation increased by 1.13±0.11 mas over the 4 epochs along a line connecting them (corresponding to a 
transverse velocity of 14.4 km s _1 ). The spatial and velocity structures of the intense group and the relative 
proper motions strongly suggest that the maser emission is associated with a protostellar jet. Comparing the 
observed LSR velocities with calculated radial velocities from a simple biconical jet model, we conclude that the 
most of the maser emission are likely to be associated with an accelerating biconical jet which has large opening 
angle of about 70° . The large opening angle of the jet traced by the masers would support the hypothesis that 
poor jet collimation is an inherent property of luminous (proto)stars. 

Key words. Stars: formation - Radio lines: ISM - ISM: jets and outflows - ISM: individual objects: IRAS 
20050+2720 MMS1 



1. Introduction 

Water maser surveys using single-dish radio telescopes to- 
ward intermediate- and low-mass young stellar objects 
(YSOs) have been extensively performed since the early 
1990s (e.g., Claussen et al. 1996). From a multi-epoch 
survey toward low-mass YSOs (bolometric luminosity, 
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Lboi < 100 Lq), Furuya et al.(2001, 2003) found that 
class objects are favored sites for the masers: the detec- 
tion rates are derived to be ^40% for class 0, while only 
~ 4% for class I. It is known that the isotropic maser lumi- 
nosity, £h 2 Oj correlates well with the bolometric luminos- 
ity of the source (Wilking et al. 1994; Furuya et al. 2001). 
It is interesting, however, to note that the presence of the 
maser emission is strictly related to that of high- velocity 
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outflowing gas in the case of high-mass YSOs (Felli, Palagi 
& Tofani 1992). |Furuya et al. (200T)| showed that the H 2 
maser luminosity in low-mass stars is more closely related 
to the luminosity of 100 AU scale radio jets rather than the 
mechanical luminosity of larger scale CO outflows. In fact, 
VLA observations have revealed that the masers tend to 
be distributed within several hundred AU of the central 
stars (e.g., Wootten 1989). Although some H2O masers 
are reported to be associated with protostellar disks (e.g., 
Fiebig et al. 1996; Torrelles et al. 1998; Seth, Greenhill, 
& Holder 2002), high resolution VLBI observations have 
demonstrated the presence of knots and shock structures 
which are reminiscent of those of ionized jets in the larger 
scale Harbig Haro objects, suggesting that in most cases 
the masers originate in shocks produced by jets from pro- 
tostars (e.g., Claussen et al. 1998; Furuya et al. 2000). 
There are a few published VLBI water maser studies of 
the jets and outflows from intermediate-mass young stel- 
lar objects (Patel et al. 2000; Seth et al. 2002). In order to 
extend our knowledge of H2O masers in intermediate- mass 
YSOs, we have conducted multi-epoch VLBA observations 
of H2O masers towards the intermediate-mass YSO IRAS 
20050+2720 (d =700 pc). 

2. IRAS 20050+2720 

IRAS 20050+2720 is surrounded by a large cluster of low- 
mass stars (Chen et al. 1997; Wilking et al. 1989) and 
has a luminosity in the IRAS bands of 388 Lq (Molinari 
et al. 1996). The IRAS source has been categorized as 
a luminous class protostar in the early compilation 
HBachiller 1996|> . Recent SCUBA imaging (Chini et al. 
2001) revealed the presence of a bright central object 
(IRAS 20050+2720 MMS1) and two associated objects 
located ~ 2' south-east. The brightest source MMS1 is 
identified as the IRAS source. Although i20 was not cat- 
egorized as class in the updated compilation (Andre, 
Ward-Thompson & Barsony 2000), pirni et al. (200ljlre- 



ported that the source MMS1 is very likely to be at the 
class stage. This is because MMS1 shows a ratio of FIR 
luminosity (Lfir) and submillimeter luminosity (L smm ) 
for A > 350 /im of ~ 164 which satisfies one of the defini- 
tions of class (Andre, Ward-Thompson & Barsony 1993; 

smm 

< 200). 

Bachiller, Fuente & Tafalla (1995; hereafter BFT95) 
found three pairs of outflow lobes emanating from the 
vicinity of the source MMS1 from the IRAM 30-m tele- 
scope CO J =2 1 observations. One of the lobe pairs is a 
highly collimated jet with extremely high velocity (EHV) 
emission whose terminal velocity exceeds ~ 60 km s _1 
with respect to the ambient cloud velocity. The presence 
of the EHV outflow suggests that the driving source of the 
EHV outflow is in its most powerful outflow phase. BFT95 
suggested that two or three independent outflows are em- 
anated from different YSOs embedded in the cloud core, 
although the driving sources have not been identified. 

H 2 maser emission in IRAS 20050+2720 region was 
first detected by Palla et al. (1991), and was subse- 



quently observed by the Arcetri group (Brand et al. 1994; 
Palumbo et al. 1994): all of the detected emission was 
seen around the cloud velocity (Vlsr = 6 km s _1 ). Using 
the Nobeyama 45-m telescope, Furuya et al. (2003) newly 
detected EHV maser emission at Vlsr = —91 km s -1 . 
The EHV emission was blueshifted with respect to the 
cloud velocity, while no high velocity emission was de- 
tected on the redshifted side. This source also showed 
weak, blueshifted, intermediate high velocity (IHV) com- 
ponents at Vlsr — —24 and —36 km s _1 in 1998 February. 
In 1999, we carried out VLA observations and found that 
all of the maser emission was located within ~ 5" (^350 
AU) from the source MMS1. Our VLA observations re- 
vealed that the EHV emission is located exactly at the 
JCMT position for source MMS1, while the low-velocity 
components around the cloud velocity are located to 4.50" 
west and 0.60" north from the EHV emission (Furuya et 
al. 2003). Suspecting multiplicity of CO outflows and the 
fact that the SCUBA beam (8.3" at 450 ^m; Chini et 
al. 2001) is larger than the separation of the two maser 
components, we made high resolution continuum images 
of the region with the OVRO mm-array. In order to inves- 
tigate the detailed structure of the masers, we carried out 
extremely high angular resolution VLBA observations. 

3. Observations 

3.1. 3 mm Continuum Emission Observations with 
OVRO mm-Array 

Aperture synthesis observations of continuum emission 
at 3 mm were carried out using the six-element Owens 
Valley Radio Observatory (OVRO) Millimeter Array 
from 2002 December to 2003 February with the H and 
E array configurations. The phase tracking center was 
a(J2000)=20 h 07 m 5"87, c5(J2000)=27°28'59"80. The field 
of view (FOV) was 65". All of the element antennas are 
equipped with SIS receivers having system noise temper- 
atures in double-sideband of 200 K toward the zenith at 
93 GHz. We tuned the 3 mm SIS receiver at the frequen- 
cies of N 2 H+ (1-0) line (93.173 GHz) for upper sideband 
and H 13 CO+ (1-0) line (86.754 GHz ) for lower side band. 
A detailed presentation of the results of the molecular line 
emission will be published elsewhere (R. S. Furuya et al. 
2005, in preparation). The Continuum Correlator was con- 
figured for both the sidebands, with a total bandwidth of 
3 GHz. We used 3C 454.3 and 3C84 as a passband cal- 
ibrator and J2025+337 as a phase and gain calibrator. 
Flux density of J2025+337 was measured by comparison 
of that of Uranus: it was stable in the range from 1.2 to 
1.5 Jy during the observation period. The overall flux un- 
certainty is about 20%. The data calibration was done 
using the originally developed software at OVRO, and 
the image construction was performed using the AIPS 
package of the NRAO. After merging the data in both 
the sidebands, we constructed continuum emission images 
with two beam weightings. Synthesized beam sizes were 

4. "04 x 3."18 with natural weighting and 2 "60 x 1 "85 



Furuya et al.: Proper Motion of H 2 Masers in IRAS 20050+2720 MMS1 
Table 1. Summary of VLB A H2O Maser Observations 



3 



MAIN Field: Low- Velocity Emission Sub-Field" : EHV Emission 



Epoch" 


Number of 


$maj X $min 


P.A. 


Sensitivity 


Number of 


Praaj X $min 


P.A. 


Sensitivity 
(mjy beam ) 




Antennas 


(mas x mas) 


(deg) 


(mjy beam -1 ) 


Antennas 


(mas x mas) 


(deg) 


I 


W d 


0.78x0.40 


-2 


3.4 


5> 


2.8x1.4 


+29 


6.2 


II 


9 e 


1.1x0.51 


+26 


3.7 


5' 


2.8x1.5 


+29 


4.7 


III 


10 d 


0.95x0.40 


-14 


2.8 


5-'- 


2.9x1.6 


+24 


4.8 


IV 


10 d 


1.1x0.39 


-14 


3.5 


5-'- 


2.9x1.6 


+27 


6.9 



(a) The 4 epochs are 1999 April 1, May 5, June 5 and July 4. (b) The following data are common for the Sub-Fields 1 and 2 
(see H4.2.2H . (c) An RMS image noise level with a velocity resolution of 0.2 km s . (d) All 10 VLBA stations, (e) Except the 
North Liberty station. (/) Fort Davis, Los Alamos, Pie Town, Kitt Peak and Owens Valley. 




20 h 07 m 06 s 05 s 4 20 h 07 m 06 s 05 s 4 

R.A.(J2000) 



Fig. 1. Continuum emission maps taken with the OVRO millimeter array toward IRAS 20050+2720 MMS1 with (a) 
natural weighting (left) and (b) uniform weighting (right). The solid contours start at +2cr level with +2cr interval 
and the dashed contours are —2a intervals started at —2a level. The image noise levels are 0.78 and 1.5 mjy beam -1 
for the panels (a) and (b), respectively. The white dot and the eastern cross indicate the positions of the low- velocity 
and EHV H2O maser emission (see H4.2.2I) . The sizes of the synthesized beams are presented in the bottom right 
corners. 



with uniform weighting. The la rms noise levels for the 
continuum emission maps were 0.78 mjy beam -1 for the 
former and 1.5 mjy beam -1 for the latter. 

3.2. H 2 Maser Observations with VLBA 

VLBI observations of the H2O maser emission in IRAS 
20050+2720 MMS1 were carried out using all 10 antennas 
of the Very Long Baseline Array (VLBA) of the NRAO 1 
on 1999 April 1, May 5, June 5, and July 4 (hereafter 
epochs I, II, III and IV, respectively). For epoch II, how- 
ever, we could not use the antenna at North Liberty. All 
of the data were obtained for 8-hour integration in each 
epoch. We used a frequency setup of the 8 MHz IF band- 
width mode with 512 channels, which provides a velocity 
resolution of 0.2 km s -1 at 22.235077 GHz. This frequency 
setup covers the range of Vlsr = —96.1 to +10.5 km s -1 : 

1 The National Radio Astronomy Observatory (NRAO) is 
operated by Associated Universities, Inc., under cooperative 
agreement with the National Science Foundation 



this velocity coverage is sufficient to detect all of the maser 
emission previously detected. 

The data were correlated at the NRAO Array 
Operation Center (Socorro, New Mexico). We adopted a 
correlator averaging time of 2.16 sec to obtain ~ 0.8" ra- 
dius FOV for the baseline of ~ 4000 km. Data calibration 
and image construction were performed using the AIPS 
package developed by the NRAO. We used two bright 
quasars 3C345 and 3C454.3 to determine delay and fringe 
rates as well as to calibrate bandpass response. In the next 
section, we present the image construction, identification 
of the maser emission and further analysis together with 
the results. 

4. Results and Analyses 

4.1. Relation between Millimeter Continuum Emission 
and Masers 

Figure ^ presents the OVRO continuum emission maps 
together with the positions of the H2O masers. There is a 
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distinct continuum emission peak toward the low-velocity 
and EHV maser emission. To plot the absolute position of 
the low-velocity masers, we adopted the results from the 
VLA measurements by Furuya et al. (2003) and used the 
position offset of the EHV emission obtained in H4.2.2I 
Clearly the millimeter continuum emission is associated 
with the low-velocity maser emission. The peak of the 
continuum emission is slightly shifted with respect to the 
maser emission. We believe that this positional shift is real 
considering the baseline accuracies of VLA and OVRO ar- 
ray and angular separations between the calibrators and 
the source. It is noteworthy that the continuum emission 
in the natural weighting map (Figure ^i) is elongated to 
the east, namely, toward the EHV emission. In the uniform 
weighting map (Figured), one may notice that there is 
a weak emission peak close to the EHV maser emission: 
this fact might suggest that there are at least two YSOs 
in this region. 

We estimate molecular hydrogen mass (Mh 2 ) of the 
core from the total flux density (F u ) of the contin- 
uum emission assuming that the emission is thermal ra- 
diation from dust grains. We used a relation Mh 2 = 
F„D 2 1 (K„B„(Td)) where k v is the mass opacity coefficient 
of the dust, Td is the dust temperature and B V (T&) is the 
Planck function. The value of k„ at 3 mm is calculated 
with the usual form of n v oc k^v^ . In order to keep con- 
sistency with the previous JCMT measurements by Chini 
et al. (2001), we used the same k of 0.003 cm 2 g -1 at 
231 GHz, /3 of 1.4, and T d of 34 K. We obtained M H2 
of 0.007Mq for F v — 25.6 mJy by integrating the emis- 
sion inside the 3<r contour. Note that the derived mass is 
a lower limit because interferometric observations do not 
receive the whole flux from a source due to the lack of 
short spatial frequency data. In fact, the resultant pro- 
jected baseline length of our OVRO observations ranged 
from 6.2 to 75 kA, which will miss 50% of the flux from 
structures extending more than 16" (0.05 pc at d =700 
pc)(see Wilner & Welch 1994). 

4.2. Maser Emission Search in the VLBA Data: Low 
and Extremely High Velocity Emission 

In the following, we present results and analyses from the 
VLBA observations of the H2O masers. 

First, we carried out fringe-frequency analysis (e.g., 
Walker 1981; Walker, Matsakis & Garcia-Barreto 1982) to 
cover a large FOV of ~ 3" which is 3 orders of magnitude 
larger than the VLBA fringe-spacings. The purpose of the 
fringe-frequency analysis was to search for maser spots 
which were not excited during the VLA observations in 
1999 February (Furuya et al. 2003). As expected from the 
VLA observations, we confirmed that the distribution of 
the low-velocity masers is sufficiently compact to perform 
standard Fourier synthesis. The EHV emission was too 
weak to be detected with the fringe- frequency analysis. 

Subsequently, we performed self-calibration using a 
strong (~ 45 Jy measured by Furuya et al. 2003) and 



point-like maser spot identified at Vlsr = +1-6 km s -1 as 
a model. In the self-calibration procedure, we solved the 
time variation of the complex gain for phase weighted by 
amplitude. 

4.2.1. Main Field: Low-Velocity Emission 

Applying the solutions from the self-calibration procedure, 
we carried out image construction toward the low- velocity 
emission (hereafter we refer to it as the MAIN Field): the 
image area has ~ 0.23" size divided into 9 fields each of 
which has 77 milli-arcsecond (mas) size and 512x512 pix- 
els with a cell size of 0.15 mas. We believe that the area 
size of ~ 0.23" was sufficient to search for maser emission 
considering the absolute position accuracy of the VLA ob- 
servations (32}. Subsequent to this coarse search, we con- 
structed a final image of 41 mas size for the area where 
the low- velocity emission was detected, with a smaller cell 
size of 0.08 mas. The image noise level per velocity chan- 
nel was typically 3 mJy beam" 1 and the synthesized beam 
size was typically 1.0x0.5 mas (Tablc^. 

In Figure El we show a total integrated intensity map 
of the MAIN Field, namely low- velocity H2O masers, ob- 
tained in the Epoch III, in which we attained the high- 
est sensitivity among the 4 epochs. The inserted panel 
presents the spectrum of the maser emission detected in 
the region. The overall distribution of the masers was sim- 
ilar for all 4 epochs: there can be seen bright maser emis- 
sion peaks at the field center (hereafter MAIN group) and 
an isolated emission peak at ~ 25 mas (corresponding to 
~ 18 AU) south-west of the MAIN group. Hereafter we 
call the latter emission as "SW feature" instead of "SW 
group" because the emission showed a point-like structure. 
The definition of "feature" will be given in H4.3I 

4.2.2. Extremely High Velocity Emission 

In addition, we searched for the EHV emission on the ba- 
sis of the snapshot VLA D-array observations in 1999 
February (Furuya et al. 2003) and 2003 January (M. 
Claussen, private communication): the former observa- 
tions (hereafter Sub-Field 1) showed that the Vlsr = — 93 
km s _1 emission is shifted by +2.33" in R.A. and —0.60" 
in Deck with respect to the low-velocity emission and the 
latter observations (hereafter Sub-Field 2) showed that 
the Vlsr = ~72 and —78 km s -1 emission is shifted by 
+4.4" in R.A. and —0.6" in Deck. Since these positions 
are outside the FOV correlated toward the 1.6 km s _1 
emission, we used UV data only from south-western 5 an- 
tennas (Fort Davis, Los Alamos, Pie Town, Kitt Peak and 
Owens Valley) which provide the maximum spatial fre- 
quency of ~ 1.1 x 10 8 A. The synthesized beam sizes were 
typically 2.8x1.5 mas and the image noise levels per ve- 
locity channel were typically 5 mJy beam -1 (Table 
Applying the solution from the self-calibration procedure 
above, we constructed images toward both the Sub-Fields. 
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Fig. 2. Total integrated intensity map of the H 2 masers 
in the MAIN Field of IRAS 20050+2720 MMS1 obtained 
at the 3rd epoch (1999 June 5) of the VLB A observations. 
Contour levels are 10, 20, 30, 50, 100, 200, 400, 800, 1600, 
3200, 6400, 12800, and 25600 times the image noise level 
of 12.3 mjy beam -1 km s _1 . The inserted panel in the 
top right corner shows the H2O maser spectrum obtained 
by integrating the emission over the presented region. The 
vertical arrow shows the ambient cloud velocity (V^ ys ) of 
Vlsr — 6.0 km s™ 1 estimated from CS J = 3-2 line ob- 
servations (BFT95). The size of the synthesized beam is 
presented in the bottom right corner. 

We detected the maser emission peaked at Vlsr = 
-93 km s" 1 in the Sub-Field 2 where the Vlsr = -72 
and —78 km s _1 emission was detected with the VLA in 
2003 January. However, we did not see any maser emission 
around Vlsr = —72 and —78 km s _1 toward the two Sub- 
Fields. The detected —93 km s _1 emission shows a point- 
like structure located at ~4184 mas east and ^709 mas 
south with respect to the Vlsr = 1.6 km s _1 emission 
peak with which we performed the self-calibration. 

Table 3. Mean LSR- Velocities (V mea n) for the Identified 
Maser Features 



Vmcan (km S L ) 



Feature 


I 


II 


III 


IV 


2 


-1.2 


-1.1 


-1.8 




4 


-0.39 


-0.30 


-0.23 


-0.1 


5 


1.6 


1.6 


1.6 


1.6 


9 


7.5 


7.0 


6.4 


6.1 


13 


-0.30 


-1.0 


-1.6 


-1.1 


SW 


3.5 


2.9 


2.9 


1.8 


EHV 


-91.2 


-91.3 


-91.2 


-91.3 



(a) Errors are typically less than 0.03 km s 



4.3. Identification of Maser Spots and Features 

To identify maser "spots" from the image, we fitted a 
two-dimensional elliptical Gaussian profile to individual 
possible spots. We adopted a detection threshold of a 
signal-to-noise ratio (S/N) larger than 10. In this way, 
we detected 30^70 maser spots at each epoch (Figure 
|3J). We estimated a relative positional error for each spot 
(fspot) to be < 0.02 mas: a relative positional error of 
a point-like source convolved with a Gaussian shaped 
beam is given by the relation of a = 0.45 # syn /[S/N] 
l|Condon 1997J) where 6 syn is the FWHM of the synthe- 
sized beam (0 syn = [6> maj x 6> min ] 1/2 )- 

Subsequent to the identification of the maser "spots" , 
we divided them into groups, i.e., spatially localized "fea- 
tures" with distinct peaks in line profiles: we grouped the 
spots in adjacent velocity channels that are distributed 
within one synthesized beam width. Each feature is con- 
sidered probably to represents a distinct clump of gas. 
In Figure |3| only Features 4, 5, 9 and 13 showed emis- 
sion over the 4 epochs among the 14 features identified 
in the MAIN-Field. Feature 2 persisted from Epoch I to 
III, but disappeared in the Epoch IV. The other 9 fea- 
tures did not persist continuously more than 3 epochs. 
Therefore, we do not consider these 9 features for our fur- 
ther analysis in proper motion measurements. For each 
feature, we calculated an intensity-weighted mean position 
of the contributing maser spots: its uncertainty is given by 
^feature = E spo t ^spotl 1/2 ■ As summarized in Tabled the 
resultant positional errors were a few xO.Ol mas for the 
6 features in the MAIN-Field and a few xO.l mas for the 
EHV feature in the Sub-Field. 

Finally, it is noteworthy that there is a velocity gradi- 
ent of approximately 9 km s _1 over 10 mas from Features 
1 to 9. 

4.4. Spectra of the Features and Their Time Variation 

For each maser feature identified, we made a spectrum by 
integrating the intensity over the corresponding region. 
Figure 0] represents spectra of the features listed in Table 
|5J Among the series of the 7 spectra, Features 2, 5 and 
EHV showed single-peaked spectra while Features 4, 9, 
13 and SW displayed multi-peaked spectra. Using these 
spectra, we evaluated intensity-weighted mean velocities 
(Vmoan; Table 01 . It is noteworthy that the former single- 
peaked features did not show any prominent trend of LSR- 
velocity change over the 4 epochs. On the other hand, the 
latter multiple-peaked features displayed small drifts. We 
discuss the velocity drifts when we assess proper motions 
of the features in the next subsection. 

4.5. Proper Motions 

Since VLBI observations generally do not provide abso- 
lute positions, we adopt the brightest spot at Vlsr = 1-6 
km s _1 in Feature 5 as a reference to examine the cross- 
epochal positional shifts of the maser features. The selec- 
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Fig. 3. Total integrated intensity maps (contours) overlaid with isovelocity maps (images) for the MAIN Field of the 
H2O masers in IRAS 20050+2720. The contour levels are the same as those in Figure [21 The velocity range shown 
here is from Vlsr = —1-5 to 8.5 km s _1 . The white dots present peak positions of the maser spots. Boxes with labels 
indicate identified maser features. 

Table 2. Relative Positions of the Identified H2O Maser Features 







Epoch I 






Epoch II 






Epoch III 






Epoch IV 




Features 


Aa a 


AS" 




Aa 


A<5 


a 


Aq 


AS 


a 


Aa 


A<5 


a 




(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


2 


-0.046 


3.928 


0.011 


0.049 


4.094 


0.014 


-0.133 


3.860 


0.012 








4 


-0.160 


1.500 


0.031 


-0.124 


1.434 


0.010 


-0.154 


1.308 


0.012 


-0.172 


1.510 


0.008 


5 


-0.133 


0.135 


0.008 


-0.034 


0.113 


0.009 


-0.053 


-0.021 


0.018 


-0.071 


-0.037 


0.019 


9 


-1.794 


-2.481 


0.014 


-1.679 


-2.436 


0.012 


-1.659 


-2.398 


0.019 


-1.684 


-2.204 


0.010 


13 


-5.439 


-1.280 


0.010 


-5.392 


-1.464 


0.016 


-5.349 


-1.600 


0.023 


-5.235 


-1.395 


0.019 


SW 


-23.11 


-12.89 


0.052 


-23.43 


-13.10 


0.04 


-23.62 


-13.50 


0.03 


-23.97 


-13.60 


0.05 


EHV 


4184.8 


-709.9 


0.29 


4184.3 


-709.9 


0.35 


4183.9 


-709.7 


0.18 


4183.6 


-709.4 


0.17 



(a) Right Ascension offset with respect to the 1.6 km s 1 spot, (b) Declination offset, (c) Position error 



tion of the 1.6 km s _1 spot is justified because this spot 
does not show any velocity drift over the 4 epochs. In ad- 
dition, the spot is enough small to be used as a position 
reference as shown in the following. The correlated flux 
(i.e., fringe amplitude) of the spot drops by a factor 2 at 



a projected baseline length of ~ 3.2 x 10 8 A which corre- 
sponds to a fringe spacing of ~ 0.06 mas. This means that 
FWHM of the maser spot would be ~ 0.06 mas assuming 
that the spot has a Gaussian shape brightness distribu- 
tion. Note that the reference spot is not a complete point 
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Fig. 4. H 2 maser spectra toward Features 2, 4, 5, 9, 13, SW and EHV in IRAS 20050+2720 MMS1. The ambient 
cloud velocity is Vlsr = 6.0 km s" 1 (BFT95). In each feature, the spectra in the left-hand panels are magnified in the 
right-hand panels. 



source compared with the relative position accuracy of the 
features (Tabled . Therefore, the size of the reference spot 
will be treated as a positional uncertainty in our analyses. 

In Figure[Sl we present an overlay of the peak positions 
of the maser spots together with boxes indicating the iden- 
tified features. The overall distribution of the maser spots 
does not change over the 4 epochs: it shows an arc-shaped 
structure. In Figure we present magnified overlays of 
the 7 features together with the intensity-weighted mean 
positions for the observing epochs (Table GJ. As expected 
from Figure most of the maser features in the MAIN 



Group do not show any prominent positional shifts. On 
the other hand, Features SW and EHV displayed distinct 
positional shifts as can be seen in Figure the member 
spots of SW sequentially appeared from NE to SW over 
the 4 epochs and those of EHV appeared from SE to NW. 
These facts strongly suggest that the observed position 
shifts are caused by real motions of the spatially localized 
masing gas. In contrast, Features 2, 4 and 5 show neither 
well-defined position shifts of member spots nor system- 
atic motions of their intensity- weighted mean positions. 
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Fig. 5. Overlay of the peak positions of the low-velocity H 2 maser spots in the MAIN-Field of IRAS 20050+2720 
MMS1. The colors present the observing epochs and the dashed-boxes indicate the maser features identified in <j4.3l 
The dashed-boxes are the same as the solid boxes in Figure [3] 



In order to identify proper motions more quantita- 
tively, we made plots of the intensity-weighted mean posi- 
tions vs. time (Figure UJ. Table 0] summarizes the results 
of our analysis: we showed the derived position displace- 
ments for the 7 features. After fifth column, we show re- 
sults only for the maser features being identified to have 
proper motions. We assess that the position shifts seen 
only in Features SW and EHV represent real motions of 
masing gas because these features displayed position shifts 
in \n\ exceeding 3cr levels of the uncertainties. In Table 
01 we present P.A. of the proper motions, transverse ve- 
locities (Vtrans) in the plane of the sky. In addition, we 
estimated 3-Dimensional (3D) velocities (T^jd) and an in- 
clination angle (i) to the plane of the sky. 

On the other hand, we conclude that the mean posi- 
tion shifts of the Features 9 and 13 were not real gas mo- 
tions, but apparent changes because their spot distribu- 



tions seem random. Let us suppose maser appearance re- 
flects the motions of shock fronts through a gas clump, as 
seems likely. Then the observations may catch the clump 
harboring Feature 9 as one shock dies out in epoch I (e.g. 
SE-NW line of maser spots in Figure [fjj) and as another 
becomes dominant in epochs II, III, and IV (e.g. NE-SW 
line of maser spots in Figure EJ). In this interpretation, 
there appears to be little proper motion between epochs 
II and III, however there does appear to be some between 
epochs III and IV. Since the status of these aspects of 
Feature 9 remain unclear, we do not include them in our 
analysis of proper motions. A dataset with more closely 
spaced observations, extending over a similar time period, 
could possibly clarify the situation. 
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Fig. 6. Overlay of the peak positions of the maser spots in Features 2, 4, 5, 9, 13, SW and EHV where the mascr 
emission was detected continuously more than 3 epochs. The colors indicate the observing epochs as in Figure [S] Each 
color-coded plus mark indicates the intensity-weighted mean position of the features at each observing epoch. The 
black cross in the Feature 5 indicates the adopted positional reference spot at Vlsr = 1.6 km s _1 . 



4.6. Results of EHV Maser Emission 



In H4.5I we presented that the observed position shift of 
the EHV maser emission is the proper motions of maser- 
ing gas. It should be noted that the proper motion of the 
EHV masers is almost directed toward the MAIN when we 
compare Figures ^ with Together with its line of sight 
velocity, we estimate that 3D motion of the masering gas 
has inclination angle of 79.4° ± 0.4° with respect to the 
plane of sky (Table |2J, suggesting that the gas motion is 
closely parallel to the line of sight. Applying the inclina- 
tion to the apparent angular separation of ~ 4.2" between 
the MAIN and EHV (Table H corresponding to ~ 2900 
AU at d = 700 pc), their 3D separation is estimated to 
be approximately 23" (~ 16,000 AU ~ 0.078 pc). Even 
if considering possible uncertainties, the real separation is 
very likely to be in the range between ~ 2,900 AU and 
~ 16,000 AU. 



We argue that the EHV masers represent a different 
protostar's activity from the Main for the following rea- 
sons. First, the separation is too large to associate the 
EHV masers with the Main, which is likely to be excited by 
the intermediate-mass protostar (§2}- In fact, Furuya et al. 
(2003) reported that all the twenty H2O maser sources de- 
tected in their VLA survey toward low- and intermediate- 
mass YSOs were associated with the central protostars 
within < 50 ~ 200 AU (see, e.g., Terebey, Vogel, & Myers 
1992). On the other hand, Homer & Churchwell (1996) 
showed that a median separation between H2O masers and 
ultra compact Hn regions is 0.1 pc for OB stars. Clearly, 
these facts support the above conclusion. Second, if masers 
are originated in the larger scale EHV CO outflows from 
the MAIN, we could have observed expanding motions. 
However, our proper motion measurements clearly showed 
that their separation has decreased. On the basis of the 
large separation and the direction of proper motions, we 
thus rule out the hypothesis that EHV masers are associ- 
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Fig. 7. Plots of intensity-weighted mean positions of the 7 features vs. time. Error bars include the positional uncer- 
tainties in Tableland the size of the reference spot at 1.6 km s -1 ( ^4.5[1 . The straight lines are the best-fit to the 
data. The origin of the time axis is 1999 April 1. 

Table 4. Position Shifts and Proper Motions for the 7 Maser Features 



Feature 


Aa a 


AS b 


M c 


P.A. d 


v trans 




V 


(mas) 


(mas) 


(mas) 


(deg) 


(km s _1 ) 


(km s -1 ) 


(deg) 


2 h 


-0.018±0.31 


-0.15±0.16 


(0.15±0.35) 










4 


0.013±0.13 


-0.017±0.011 


(0.021±0.14) 










5 


-0.17 ±0.034 


0.042±0.03 


(0.18±0.076) 










9 


0.28±0.058 


0.098±0.043 


(0.30±0.094) 










13 


-0.13 ±0.017 


0.19±0.062 


(0.23±0.06) 










SW 


-0.75±0.078 


-0.85±0.085 


1.13±0.11 


-131.4 


14.4±1.4 


14.4±1.4 


4.7±3.5 


EHV 


-0.50±0.12 


±1.26±0.085 


1.36±0.16 


-68.2 


17.9±5.2 


-91.3±29.5 


79.4±0.2 



(a) and (b) Total position shifts of the maser features over the observing period, namely per 95 days, along R.A. and Deck 
directions, respectively, derived from the least-square fitting in Figure □ (c) = ^(Aq) 2 ± (AS) 2 . For the maser features 
where we could not detect well-defined proper motions, apparent position shifts are given in the parenthesis, (d) Position Angle 
of the proper motion vector in the plane of the sky. (e) Transverse velocity on the plane of the sky converted from (/) 
3-Dimensional velocity obtained from V-jd = Varans ± (Kncan — V TC t) 2 where Kioan is in Tableland V rc f is the LSR-velocity of 
the reference spot (1.6 km s -1 ). Negative and positive represent approaching and receding motions, respectively, (g) Inclination 
angle of the V3D vector with respect to the plane of the sky, obtained from i = tan -1 (| V mC aii — Kef l/Krans). (h) All the parameters 
for Feature 2 are derived from the data taken in the first 3 epochs. 



ated with the jet emanated from the MAIN. We suggest 
that the EHV masers must be associated with another 
member source in the cluster (Chen et al. 1997; Wilking 
et al. 1989). 

In order to search for the possible exciting source of the 
EHV maser emission, we analyzed 4.86 GHz radio con- 
tinuum emission data from the VLA Archive Database. 
We, however, could not detect a compact radio contin- 
uum emission with a 3tr upper limit of 0.14 mJy beam -1 
(flmaj x min ~ 5.97" x 5.15"). On the basis of the lack 
of a bright compact continuum source and the presence 
of the EHV maser emission, we speculate that the driving 
source of the EHV emission could be an extremely young 
protostar. 



5. Discussion 

Evidence from AU-scale VLBI H2O maser observations 
suggests that the masers in star forming regions are most 
likely to be excited in the interaction zone between a jet 
and ambient cloud material; or in interaction between a 
wide angle flow and the surface of a protostellar disk. 

In this section, we discuss the origin of H2O masers 
associated with the millimeter continuum source MMS1 
in terms of the jet scenario, supported by the following 
evidence. We do not further discuss the EHV maser emis- 
sion which we have shown to be associated with other 
YSO activity ( i)4.6fl . We first note that the spatial and 
velocity structures of Features 2 to 12 in the MAIN group 
(Figure |2| convincingly demonstrate that they are ex- 
cited in outflowing gas associated with a protostellar jet. 
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+x 

Fig. 8. Sketch of the biconical jet model. The z and x axes 
correspond to the line of sight and the jet axis projected 
on the plane of the sky, respectively. An observer lies at 
z = — oo, and the x — y plane corresponds to the plane of 
the sky. 
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Fig. 9. Guide map showing the 3 jet axes that have been 
taken as cutting lines in the position-velocity diagrams 
(see text; Figure HU)l . To avoid confusion, we present only 
an outline of each single lobe of the jet. The position angles 
of the 3 jet axes in the upper panel are 61° for the axis 
with the labels of Al and A2, 101° for Bl and B2, and 
28° for Cl and C2. 



The maser velocity gradient parallels a line connecting 
MAIN and SW masers. It also parallels the NE-SW pair 
of the Intermediate High Velocity (IHV) CO outflow lobes 
(P.A.~ -150°; BFT95), although on a smaller scale than 
the 10" (corresponding to 0.034 pc) scale of the CO lobes. 
Secondly, the velocity sense of the two flow signatures 
agrees: the blueshifted masers lie to the NE side and the 



redshifted masers lie to the SW, the same as found in the 
velocity structure of the IHV CO outflow lobe pair. Third, 
the relative proper motion between the MAIN and SW 
shows expansion, and its direction (P.A.= —131°; Table 
0} is almost parallel to the flow line. These results strongly 
suggest that the maser jet channels along the direction of 
the line connecting MAIN and SW. It would be difficult 
to reconcile these motions with interaction of a flow with 
a protostellar disk. 

We assume that the exciting source of the jet is located 
very near the position of the reference 1.6 km s^ 1 spot for 
the following reasons. First, Feature 5, which hosts the 
1.6 km s _1 spot, shows a single-peaked spectrum over all 
4 epochs (Figure |3J , and it did not show a velocity drift 
(Table 0). Second, the intensity weighted mean position 
of Feature 5 and position of the 1.6 km s _1 spot showed a 
positional coincidence within 0.11 mas over the 4 epochs 
(see FigureEJ, which is consistent with the result that the 
reference spot would have a FWHM~ 0.06 mas derived 
from the fringe amplitude analysis ( H4.5[l . Last, the inten- 
sity weighted mean positions lie near to a line connecting 
the MAIN and SW masers, suggesting that the 1.6 km s _1 
spot is located at the expansion center. 
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Fig. 10. Comparison of the radial velocities from the 
biconical jet model with the observed LSR-velocities in 
the position- velocity (PV) diagram along the A1-A2 axis. 
Here we plotted all of the maser spots detected in the 
4 epoch observations. The vertical and horizontal dashed 
lines, respectively, represent the positions of the reference 
1.6 km s _1 spot. Expected blue- and redshifted radial ve- 
locity regions from the model are hatched with blue and 
red, respectively (see text). 

Now we try to shed light on the nature of the masers 
using a simple jet model. We assume that the H 2 maser 
emission is excited in the material at the interface be- 
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tween the protostellar jet and the ambient gas, namely, at 
the surface of the cone whose opening angle is 29. Such a 
model has successfully explained the distribution of H2O 
masers in the high-mass (proto)star IRAS 20126+4104 
(Moscadelli, Cesaroni, & Rioja 2000) Figure |H1 schemat- 
ically shows the model: the assumed protostar lies at the 
apex of the cone. The axis of the cone is inclined by an 
angle of <f> with respect to the plane of the sky. We define 
a coordinate system whose z and x axes are, respectively, 
parallel to the line of sight and the projection of the jet 
axis on the sky. To calculate the jet velocities seen by 
an observer (v z ) at z ~ —00, we consider that the gas 
moves along straight lines passing through the apex into 
two opposite directions. Given a power-law velocity profile 
of v{r) — vo(^) a , the jet velocity along the line of sight 
can be written as 



v z (x) = v 



R cos( 



(1) 



Here we adopt vq of 14.4 km s _1 from the 3D-velocity 
of the SW maser (Tabled}, and its distance of 21 AU from 
the 1.6 km s _1 spot for R. We believe that this assump- 
tion is reasonable when we considered other results from 
VLBA proper motions studies in low- and intermediate- 
mass YSOs (IRAS 05413-0104: 64±22 km s^ 1 at a dis- 
tance of 40 AU from the expansion center [Claussen et 
al. 1998], S106 FIR: 25-40 km s" 1 at 25 AU [Furuya et 
al. 2000], IRAS 21391+5802: £ 42 km s" 1 at - 20 AU 
[Patel et al. 2000], NGC 2071 IRS3: 22-42 kr 1 at 
260 AU [Seth et al. 2002]). We thus have the following 
four free parameters: P.A. of the projected jet axis to the 
plane of the sky (i.e., P.A. of x-axis), </>, 9, and a. By 
definition, we can give constraints of < 9 < 7r/2 and 
< |0| < 7r/2. To apply such model, we selected 3 pos- 
sible axes of P.A.= 61°, 101°, and 28° which we refer to 
as, respectively, A1-A2, B1-B2, and C1-C2 (see Figure 
IHJ). Note that 61° is the P.A. of the line connecting the 
reference spot and the SW maser, and that B1-B2 and 
C1-C2 are parallel to the two pairs of CO outflow lobes 
(see Figure 3 of BFT95). As for the power-law indices 
of the jet velocity, we considered representative values of 
a = 0, +1, and —1 which characterize constant velocity, 
accelerating, and decelerating jets, respectively. 

We compared the observed and calculated velocities in 
position-velocity (PV) diagrams: the black dots in Figure 
ITUl represent the PV distribution of the masers along 
the A1-A2 axis. Here we took the LSR- velocity of —1.6 
km s _1 as a mean velocity. Figure ITUl clearly tells us that 
the systemic velocity of Vlsr = 6 km s _1 (see the inserted 
panel in Figure^) which was derived from the 0.01 pc scale 
molecular cloud (BFT95) is not valid for the 10 AU scale 
maser emitting region. We started searching for the best- 
fit parameters from a. We found that only accelerating jet 
can explain the observed velocity structure whereas both 
constant velocity and decelerating jets cannot. We thus 
take a = +1, leaving us with two free parameters — <f> 
and 9. We obtained the best-fit parameters of <p = —39°, 



and 9 — 36°. The blue- and red hatched regions in the 
1st and 3rd quadrants of Figure^] respectively, show the 
expected LSR-velocity ranges for the blue- and redshiftcd 
components. Although the observed maser emission is not 
excited all over the expected regions, a large opening an- 
gle such as 29 p> 70° is required to explain the velocity 
structure in the 1st and 3rd quadrants. This geometry 
suggests that the line of the sight would match the cone 
surface (that is, \9\ + \<f>\ ~ -|). However, the remaining 
maser spots in the 4th quadrant, which are Features 13 
and 14 (see Figure |3J), do not reconcile with the jet model 
prediction. 

Since only the accelerating jet explains the velocity 
structure, we extended our analysis to the remaining two 
axes of B1-B2 and C1-C2, keeping a = +1. The best-fit 
parameters for the B1-B2, and C1-C2 cases are essen- 
tially the same as those obtained from the A1-A2 in the 
following three senses: 

— (i) Requiring a large opening angle (29 = 12° for Bl- 
B2 and 74° for C1-C2), in other words, the P.A. would 
have an uncertainty of 35°. 

— (ii) The cone is most likely to have a geometry such 
that the surface where masers are excited is nearly 
parallel to the line of sight (|0| + |0| = 78° for B1-B2, 
and 80° for C1-C2), namely, the SW masers lie at the 
either of the two edges of the cone. 

— (iii) The LSR-velocities of the spots in Features 13 and 
14 cannot be simultaneously reproduced with other 
spots. 

In addition, we performed further analysis by changing 
the P.A. through 5° increments, realizing the following 
points which might jeopardize the above conclusion: 

— (a) There is a 4 order of magnitude difference in the 
spatial scales between our VLBA images and CO (2-1) 
maps. 

— (b) There is no clear association of CO outflows with 
known driving source(s) ( H4.ll BFT95; Chini et al. 
2001). 

— (c) The maser emission does not trace the whole of the 
outflowing gas. 

Nonetheless, the analysis with 5° P. A. steps showed 
that the above results from (i) to (iii) are valid, and that 
no other range of P.A. than from 25° to 110° can be taken. 

Of the available observations of this region, our high 
resolution data has the best opportunity to discern indi- 
vidual sources of outflow. We stress that a single jet model 
cannot explain the PV structure of the Features 13 and 
14 which were located at the most western portion of the 
chain of the MAIN maser features, but showed a clear 
velocity gap with respect to the coherent velocity struc- 
ture from the Features 3 to 12 (see Figure Can they 
be used to pinpoint additional sources of outflow, perhaps 
associated with the other jet-like CO outflows from the 
region? We examined the velocity structure of the discor- 
dant spots in terms of a "multiple jet scenario" by apply- 
ing the above "single jet model" again on "residual PV 
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diagrams" . We subtracted maser spots in the hatched re- 
gions in Figure 221 In addition, we subtracted all the SW 
spots, assuming that such spots are associated with the 
"main" jet. Hence, the "residual spots" consist of all the 
spots in Features 13 and 14, and most of the spots in 
Features 1-3. "Residual PV diagrams" were made along 
position angles incremented by 30°, including all the CO 
outflow position angles. Given the same driving source 
as the "main" jet, one of the second jet lobes must lie to 
the NW, which might be characterized with a P.A.~ 150°, 
26 ~ 130°, and <f> ~ 60°, but with no symmetrically placed 
lobe. We are unable to convincingly associate the resid- 
ual spots with any particular flow geometry or additional 
discrete source. We speculate that the residual spots are 
being excited in some shocked regions associated with the 
"main" jet. 

In conclusion, we summarize our jet model analysis by 
combining the results from i) to iii) with those from the 
"residual PV" analysis: 

1. The majority of the H2O masers are most likely to be 
associated with an accelerating biconical jet whose pro- 
jected axis has a P. A. of 60° ±35°, and whose opening 
angle is larger than about 70° . 

2. The cone surface enclosing the jet almost matches the 
line of sight. 

3. We could not clearly explain the origin of remaining 
maser spots, and speculate that they are being excited 
in more complex jet-cloud interactions. 

We suggest that MMS1 is probably a single star and 
is the source of the flow exciting the masers lying near 
it. If there has been no precession(s) of the maser jet and 
CO outflow(s), our results indicate that the maser jet is 
most closely associated with the NE-SW pair of the IHV 
CO lobes because their axes are almost parallel. The large 
opening angle of ~ 70° at a distance of 21 AU from the 
assumed expansion center is similar to that measured in 
the intermediate-mass YSO of IRAS 21391+5802 by Patel 
et al. (2000)-these authors reported an opening angle of 
~ 110° at the 150 AU point from the star. Interestingly 
these results are also consistent with those from high- 
mass (proto)stars (Shepherd, Claussen, & Kurtz 2001, 
and references therein) rather than low-mass protostars 
where highly collimated CO outflows are generally seen 
(Richer et al. 2000, and references therein). We posit that 
higher luminosity YSOs tend to have poorly collimated 
jets, and suggest that a similar jet morphology obtains for 
intermediate-mass stars. The conclusion that the masers 
accelerate agrees with observations of other intermediate- 
mass YSOs such as S106 FIR (Furuya et al. 1999) and 
IRAS 21391+5802 (Patel et al. 2000). In addition, this 
conclusion does not disagree with that from high-mass 
star forming region W49N where Gwinn, Moran, & Reid 
(1992) reported that expansion velocity of the masering 
gas has a constant velocity of ~ 18 km s _1 up to 0.1 pc 
in a distance from the center, beyond which their velocity 
increases to more than 200 km s _1 . 



Given the complexity of the region seen in the CO 
maps (BFT95) and the near-IR image (Chen et al. 1997), 
we believe that sub-arcsecond resolution imaging of the 
CO outflow and continuum emissions with radio interfer- 
ometers will associate jets and outflows with their driv- 
ing sources. These interferometric observations will fill 
the spatial resolution gap between our AU scale VLBI 
view of the masers and the 0.01 pc scale single-dish tele- 
scope view of the CO outflows. As we have mentioned, 
i20 has displayed H2O maser emission at Vlsr = —72 
and —78 km s _1 since 2003 January, which were not seen 
during our observations f H4.2.2l) . Subsequent monitoring 
observations using the Green Bank 100-m telescope (A. 
Wootten, private communication) showed that these emis- 
sions have flared up to ~20 Jy in 2004 October. Together 
with sub-arcsecond interferometric observations, further 
VLB A H2O maser study will help to assess the nature 
of this multiple jet-outflow system, which harbors some 
of the highest velocity outflowing gas in any star forming 
region known to date. 

6. Summary 

We have performed a monthly 4-epoch VLBA observa- 
tions of the H2O masers in the intermediate- mass pro- 
tostar IRAS 20050+2720 MMS1 together with aperture 
synthesis observations of 3 mm continuum emission with 
the OVRO array. The main results of this study are sum- 
marized as follows. 

1. From the VLBA images taken with all the 10 antennas, 
we found the two groups of the low- velocity H2O maser 
spots toward the bright millimeter continuum emission 
peak. One group (the MAIN group) showed intense 
emission around the cloud velocity, the other (the SW 
feature) was located at a projected distance of 18.2 AU 
south-west of the MAIN. The OVRO 3-mm images 
clearly showed that a bulk of millimeter continuum 
emission (MMS1) is associated with the low- velocity 
masers. 

2. Using only the south-western antennas of the VLBA, 
we have succeeded in detecting the EHV maser emis- 
sion blueshifted by 99 km s _1 to the cloud velocity. The 
EHV emission is located at 4400 mas in east and 709 
mas south of the low-velocity emission. Considering 
the large 3D separation between the MAIN and EHV 
features, and the proper motion of the EHV toward the 
MAIN, we concluded that the EHV emission is not 
associated with the low-velocity emission. A distinct 
millimeter continuum source appears to be associated 
with the EHV masers, which we speculate is the driv- 
ing source of the EHV emission and which could be an 
extremely young protostar. 

3. The overall structure of the masers in the MAIN group 
did not change over the 4 epochs. On the other hand, 
we found that the projected separation between the 
MAIN group and the SW feature increased by 1.13 
mas over the 4 epochs, which corresponds to a trans- 
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verse velocity of 14.4 km s _1 . This increment of the 
separation indicates proper motion(s) of spatially lo- 
calized masing gas in the MAIN group and/or the SW 
Features. 

4. From the analysis of the velocity field of the masers, we 
conclude that the majority of the H 2 masers in i20 
MMS1 are likely to be associated with an accelerat- 
ing biconical jet whose opening angle is approximately 
70° at a distance of 21 AU from the central star. The 
presence of such accelerating jet indicates that the cen- 
tral protostar is driving the powerful jet. Moreover, the 
obtained large opening angle of the jet would support 
the hypothesis that poor jet collimation is an inherent 
property of luminous YSOs. 
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